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We present a new class of fluctuation relations, to which we will refer as Fluctuation Relations for Current 
Components (FRCCs). FRCCs can be used to estimate system parameters when complete information about 
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sense that they do not depend on some basic types of electron interactions and some quantum coherence effects. 
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I. INTRODUCTION 

Fluctuation theorems 0] El are fundamental results in the 
nonequilibrium statistical mechanics. Their discovery led to 
optimism that they might serve as universal laws that had long 
been missing from the study of nonequilibrium systems. 

Fluctuation relations for currents [3| play a special role in 
the physics of non-equilibrium transport. Measurements of 
statistics of heat production generally involve tracing com- 
plete system's stochastic trajectories. This is often an ex- 
perimentally challenging problem. Measurements of currents 
are relatively simple. Formally, a current can be found just 
from the knowledge about an initial and a final state of a 
system, and no separation of the measured system from the 
heat bath is required. Moreover, experimental studies of fluc- 
tuation theorems typically need large statistics of events be- 
cause any Gaussian distribution of any variable trivially satis- 
fies a standard fluctuation relation. Hence, nontrivial fluctua- 
tion relations can be identified only if statistics of events are 
obtained beyond the domain of Gaussian fluctuations. Non- 
Gaussian statistics in nanoscale conductors is due to the shot 
noise of electrons at non-equilibrium conditions [4|. Recently, 
it became possible to experimentally study such non-Gaussian 
fluctuations in mesoscopic electric circuits 0. 

In application to an electric circuit with two lead contacts, 
the Fluctuation Theorem predicts that the probability distri- 
bution P[q] of observing a charge q passed between two lead 
contacts with a voltage difference V satisfies the law (3J: 

P[q]/P[-q] = e F \ (1) 

where F = V/k^T. Here T is the temperature, and fee is 
the Boltzmann's constant. Eq. ([TJ is expected to be universal, 
i.e. it should be valid independently of the type of electron 
interactions in a conductor. Surprisingly, recent experimental 
work |6| has shown that the law ([T| can fail in an electric cir- 
cuit, but could be salvaged under the experimental conditions 
of (5) if the parameter F is suitably renormalized by a factor 
~ 10" 1 . 

The need to modify ([T]l was qualitatively explained in J6) 
by presence of a feedback between measured and measuring 
circuits. For example, in experiment O, the nanoscale cir- 
cuit was connected to a read-out circuit made of an additional 
tunnel junction, which was coupled to its own leads. When 



interactions between measured and measuring currents are in- 
volved, the fluctuation theorem is applicable only to the total 
system that includes both the studied circuit and the measuring 
one but this does not imply ([T} for a single current component 
any more. 

This argument only partly explains the experimental result 
in [6 1 because it makes unclear why, after considerable renor- 
malization of parameters, an individual current through a spe- 
cific lead contact again satisfies ([lj. For example, according 
to Crook's equality, the heat W, dissipated by a complete sys- 
tem, satisfies the relation 

P[W]/P[-W] = exp (W/k B T) . (2) 

We can introduce the vectors, q and V, whose components 
are, respectively, the numbers of electrons passed through in- 
dividual lead contacts, and individual voltages at correspond- 
ing leads. Numbers of vector components is equal to the num- 
ber of independent lead contacts in the circuit. The dissipated 
heat can then be expressed as W = V-q, which corresponds 
to the standard fluctuation relation for multicomponent cur- 
rent J3): 

P[q]/P[-q] = cxp (V • q/fc B T) . (3) 

Apparently, a specific current component is not propor- 
tional to the dissipation function W because the latter depends 
on all current components. Consequently, although the full 
vector of currents satisfies Eq.([3j there seems to be no rea- 
son why an individual current through a specific lead contact 
should satisfy a fluctuation relation ([T} in a multicomponent 
system. 

The appearance of fluctuation relations ([T]l for specific cur- 
rent components at renormalized values of parameters re- 
mains poorly understood. Recent studies [7| showed that 
Fluctuation Relations for Current Components (FRCC)s, i.e. 
for currents through specific links of a circuit, can appear in 
some limits of a model that corresponds to a 4-state Markov 
chain kinetics. However, generalizations of this result have 
been unknown. Another possibility to explain experimentally 
observed FRCC was based on separation of time scales [6|, 
i.e. if one current component is macroscopic in comparison to 
another one, the former can be considered as part of environ- 
ment so that one can introduce effective dissipation function 
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for the second, microscopic, current component. Such an ex- 
planation, however, can be justified only in very specific limits 
and relatively simple circuit geometries. 

In this article, we prove that there are, in fact, general con- 
ditions under which FRCCs should hold for currents through 
certain links in complex mesoscopic electric circuits coupled 
to multiple lead contacts. Unlike standard fluctuation theo- 
rems, there is no direct relation between FRCCs and the dis- 
sipated heat in a system. Instead, we will show that FRCCs 
follow from the observation that statistics of particle currents 
depends only on probabilities of single particle geometric tra- 
jectories while the information about time moments, at which 
particles make transitions along such trajectories, is irrelevant. 
Then, there can be purely topological constraints on contribu- 
tions of geometric trajectories to currents through some links 
of a circuit. For example, if a link does not belong to any 
loop of a graph, and if transitions through this link in op- 
posite directions are counted with opposite signs, then any 
geometric trajectory can make only ±1 or valued currents 
through such a link. We will show that this restricts statistics 
of currents through this link to satisfy Eq. ([TJ even when there 
are no detailed balance constraints on kinetic rates. Another 
interesting topological constraint appears when there is only 
one reservoir that supplies/absorbs particles to/from the sys- 
tem. Then all single particle geometric trajectories have to be 
cyclic, and one can make a correspondence between currents 
through links and independent cycles that a trajectory makes. 
Probabilities of independent cycles are known to satisfy re- 
lations of type ([TJ [3 1, which eventually results in FRCCs 
in such a circuit. We will also show that FRCCs are robust 
against adding important many-body interactions because the 
latter influence timing but do not change relative probabilities 
of geometric trajectories. 




FIG. 1. Chaotic cavity coupled to N electron reservoirs at different 
potentials and temperatures. 

The structure of our article is as follows. In Section 2, we 
will illustrate FRCCs in the model of a chaotic cavity, shown 
in Fig. [JJ which frequently appears in studies of counting 
statistics JSJ . In Section 3, we will increase the complexity 
of the electric circuit geometry to demonstrate the ubiquity of 
FRCCs in mesoscopic electronics. Sections 4 and 5 play a 
supplementary role. They explore both presence and absence 
of FRCCs in specific models, which statistics of currents can 



be studied in detail. They demonstrate that exactly solvable 
models produce results in agreement with more general the- 
ory developed in Sections 2 and 3. 

H. THE MODEL OF CHAOTIC CAVITY COUPLED TO N 
LEAD CONTACTS 

The model assumes that N large reservoirs at different po- 
tentials exchange electrons with a mesoscopic conducting re- 
gion (cavity), in which the electron motion is randomized 
and is influenced by exclusion interactions due to the Pauli 
principle and Coulomb interactions J8). The model also 
assumes that time-scales for self-averaging are much faster 
than the time-scale at which number of electrons in a cav- 
ity changes with time so that we can treat interactions with a 
mean field approximation. Electrons enter the cavity through 
the lead contacts with kinetic rates, k] n (Q) = hiQifx{Q), 
and leave the cavity with rates k° ut (Q) — gif 2 (Q), where 
i = 1, . . . , N. The parameter hi depends on thermal equi- 
librium conditions in the i-th reservoir, ^ is a strength of the 
cavity coupling to the i-th contact, and Q is the instantaneous 
number of electrons in the cavity. The functions fx (Q) and 
/ 2 (Q) describe the effect of many -body interactions on ki- 
netic rates. 

We assume that all parameters, as well as functions fx and 
/2, may depend on constant temperatures and potentials at 
reservoirs. As an example, first consider that the cavity is 
small (the quantum dot limit), so that it has only a single elec- 
tron level at energy E. Coulomb interactions forbid to have 
more than one electron inside. The kinetic rate of escape 
of electron from this cavity into the zth contact can be esti- 
mated by the golden rule: fc° ut = 2TrQp i (E)\T i (E)\ 2 (l - 
1/[1 + e (s-v 1 )/(fe B T)]^ where Q = 1 or Q = depending 
on the presence or absence of an electron inside the quan- 
tum dot. Ti(E) is the element of the scattering matrix be- 
tween the state inside the quantum dot and a state in the 
reservoir at the same energy E, Pi{E) is the density of en- 
ergy levels inside the ith contact near energy E, and the 
factor 1 - 1/[1 + e ( - E ~ Vi)/(feBT:i ] is due to the Pauli princi- 
ple that forbids transitions into the filled states of the reser- 
voir. Similarly, the kinetic rate of transitions from the ith 
reservoir into the quantum dot can be estimated as k° ut = 
2n{l-Q)p l (E)\T l (E)\ 2 (l/[l+e iE - Vi '>^ kBT, >}). Wecanthen 
identify 9i = 2TT Pl (E)\T t (E)\ 2 (l - 1/[1 + efa-vOAfeT)^ 
hi = e-^-^y^^, fx = (1 - Q), and f 2 = Q. An- 
other limit of the cavity model corresponds to a large cav- 
ity with kinetic rates induced by thermal over-barrier tran- 
sitions and number of states inside the cavity much larger 
than the number of electrons. In such a classical limit, ki- 
netic rates are given by Arrhenius form kf 1 ~ e ( v i- w i)/ k ^ T ( 
and k° ut ~ e (M(Q)-w,)/fc B T where w . is size of the bar . 

rier that separates the cavity from the ith reservoir and /i(Q) 
is the chemical potential of electrons inside the cavity. Then 
we can identify g 4 = e - w ^ k * T , hi = e v ^ ksT , A = 1 and 
f-2 = e M ^)/ feBT . Note that, in both cases, at constant temper- 
ature, parameters hi satisfy the relation 

hi/hi = exp{(U, - V 3 )/k B T}, (4) 
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which guarantees the presence of the standard fluctuation re- 
lation This restriction, however, will not play any role in 
our following discussion. 
Let 

Z( x ) = e s <x) =£p(g) e *-x (5 ) 

be the generating function of currents through all leads. Here 
components of the TV-vector, q, are number of electrons that 
pass through corresponding contacts during the observation 
time. In the rest of the article, all introduced vectors will cor- 
respond to sets of elements indexed by the indices of the cor- 
responding reservoirs, e.g. x = (Xi> • ■ ■ Xn) is the vector of 
counting parameters. S{x) is called the cumulant generating 
function because its knowledge corresponds to knowledge of 
all cumulants of the current distribution, e.g. 

(&) = (dS/dxi) x =o, var( ft ) = (d 2 S/dx 2 ) x =o, etc. 

(6) 

We will detect FRCCs by setting all Xu where i 7^ k for some 
k G 1, . . . , N, to zero and observing the symmetry, 

S(0,...,0, X k,0,-,0)=S(0,...,0,-Xk+F k ,0,...,0), (7) 

where F k , for a given k, is a constant parameter. By applying 
the inverse Legendre transform, one can verify that (|7]i leads 
to the FRCC, 

P(q k )/P{-q k ) = exp(g fc Ffc). (8) 

Let us calculate the probability, of that the nearest en- 
tering electron will have a geometric trajectory that enters the 
cavity via the contact i and leaves the cavity via the contact j. 
Kinetic rates for this particle depend on Q, which may change 
with time arbitrarily, however, the ratio of either two in-going 
or two outgoing kinetic rates for the given particle remains 
constant, which means that the probability of the geometric 
trajectory for this particle can be found explicitly: 

<L k =l h k9k 2^r=l 9r 

Let P(n) be the probability that during a large observation 
time exactly n electrons enter the cavity via any node. The 
probability of a geometric trajectory in (|9]l is independent of 
n and of other particle's trajectories so the generating function 
Z(x) is given by 

00 / n \ ™ 

n=0 \*,3'=1 / 

Although it is impossible to derive an explicit expression for 
P(n), one can check that the symmetry ^ is the symmetry 
of each term in (jTOJ with 

F k = In ( ) , k = l,...,N, (11) 

\hkl^i^ k 9ij 



which proves the FRCC for the cavity model. 

Note that parameters F depend not only on lead charac- 
teristics h, which can be externally controlled, but also on 
coupling parameters g. The latter may nontrivially depend on 
voltages due to appearance of screening charges in the vicin- 
ity of components of the nanoscale circuit at nonequilibrium 
conditions 1 9 1 . Estimation of parameters g from knowledge of 
lowest current cumulants may be difficult because functions 
fi(Q) and fa{Q) influence lowest current cumulants. Sur- 
prisingly, the FRCCs do not depend on interactions encoded 
in functions f±(Q) and f%{Q) at all. Hence measurements 
of FRCCs can provide us with a unique approach to measure 
the vector g in a nonequilibrium regime irrespective of many- 
body interactions inside the cavity. 

We also note that our derivation of the generating function 
is valid only for very large observation time, so that we could 
disregard the trajectories of electrons that entered but did not 
leave the cavity. Thus, FRCCs must be understood only in the 
sense of the dominating exponent (also known as the Large 
Deviation Function) of the probability distribution of a cur- 
rent. 



III. STOCHASTIC TRANSPORT ON NETWORKS 

To explore general principles that lead to FRCCs, we con- 
sider a generalization of the cavity model to a network of 
chaotic cavities coupled to lead contacts and to each other, 
e.g. as shown in Fig. [2] Electrons enter cavities through the 
reservoirs (with the rate hiQi for the i-th cavity). If two cavi- 
ties (graph nodes) i and j are connected, each electron in the 
node i can make a transition to node j with rate 5^ ; generally, 
7^ gji. Eventually, electrons leave the network through 
one of the contacts. Physics of incoherent effects due to the 
shot noise and thermal Johnson-Nyquist noise in continuous 
conductors can be obtained from the continuous limit of such 
network models [ 10 1 . 

For simplicity, we assume that all electrons are non- 
interacting, i.e. ff = 1 and j\ = Qi for all cavities, although 
generalizations to local interactions (/^(Qi) 7^ Qi) are pos- 
sible because such interactions do not change relative prob- 
abilities of geometric trajectories iPTD . It was shown in ifTTI 
that finding statistics of currents in such a model reduces to 
solving a finite set of coupled linear algebraic equations. Al- 
though explicit solutions are bulky and not illuminating, we 
used them to check the presence/absence of an FRCC for any 
link of a network numerically, as we explain in Section 5. As 
expected, we did not find an FRCC when there was no previ- 
ously mentioned topological constraints on geometric trajec- 
tories. However, FRCCs were obtained in two wide classes of 
links with such constraints. 



A. Links that do not belong to any loop of the network 

These links represent lead contacts and also internal links 
of the network that on removal break the network into disjoint 
components. We marked such links by green color in Fig. [2] 
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Suppose that we can trace the geometric trajectory of a single 
electron. When electron enters the network it makes a single 
transition through a lead contact in positive direction and then 
makes one transition through one of the leads in negative di- 
rection when it leaves the network. Let p k j be the probability 
that the electron that enters through the contact k leaves the 
network through the contact j. If an electron enters through 
the contact fc, then the moment generating function of cur- 
rents through the contact k that this electron produces during 
its life-time is given by 



Zk{Xk) = Pkk + ^2pkje Xk . 

&k 



(12) 



The generating function of currents in the contact k, which is 
produced by an electron that enters from contact j, j ^ k, is 
given by 



Z j(Xk) = X! /'.•,•' +P]k£ 

J'^k 



Xk 



(13) 



hi 



En 2 





FIG. 2. Circuit of coupled chaotic cavities and electron reservoirs. 
Green color marks links that carry currents that satisfy FRCCs. Red 
color marks links with currents that violate FRCCs. 
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(t) be the probability that during time t exactly 
qi, . . . , qN particles enter the network through the contacts 
1, . . . , N. Since particles enter independently, this distribu- 
tion is, in fact, Poisson and its generating function is given 
by 

Z P {s) = ]T P q e*« = cxp h i9it(e Si - 1) J • (14) 

q \i=l J 

Since electrons do not interact, probabilities P qi .... : q N (t) and 
Pkj are not correlated. The statistics of currents through con- 
tact k during the whole process is then given by 



Z(Xk) 



la. 



V9l 



N > 



(15) 



where we used the fact that the generating function of a sum 
of independent processes is the product of generating func- 
tions of individual processes. The latter expression shows 



that Z(xk) coincides with Zp(s) up to a change of variables 
e Si — > Zi{xk) for any i. This fact and the explicit expressions 
for Zp(s) and Zi(xk) lead us to the conclusion that the total 
generating function of currents depends on \k only through 
the combination, hkgkJ2 j7 t k Pkje Xk + )J ,,,//,(/,/*,/ ( x '. 
which has the symmetry under the change of variables, \ k — > 
-Xk + F k, where 



F k = In 



Ej^fc h i9jVjk 

hkgk T,j^ k Pkj 



(16) 



The case with a link that connects two otherwise disjoint 
components, a and 7, is proved similarly: Let p a be the 
probability that a given particle that enters the network in 
the components a leaves the network through some contact 
in another graph component 7. The generating function of 
currents through the link that connects such components is 
Z a (x) = (1 ~ Pa) + Pa& x - The rest of the proof is the same 
as for the currents through lead contacts, where instead of con- 
tact indices k and j we write component indices a and 7. 



B. Quantum coherence among trajectories 

So far we assumed the lack of quantum interference among 
particle trajectories. The arguments leading to |7| and (16 1 
for a particle's motion on a network, however, do not refer 
to any classical or thermodynamic reason. For example, one 
can imagine that the distribution of electrons in leads is not at 
equilibrium and that probabilities pij are influenced by quan- 
tum interference of different trajectories that connect nodes i 
and j. We just should assume that (a) events of particle's es- 
capes into reservoirs destroy coherence, (b) particles enter the 
network according to the Poisson statistics, and (c) transition 
probabilities p^ are constant. Such conditions are realized in 
the quantum regime when the single electron scattering am- 



plitude, Sij(E), between any pair of different reservoirs i and 
j in a channel with any energy E, is small, i.e. \sij(E)\ <C 1. 
In this limit, we can disregard simultaneous multi-electron 
scattering processes, although purely quantum effects, such 
as quantum interference of trajectories, can still influence sy. 

The generating function for N terminals and non- 
interacting fermions was derived by Levitov and Lesovik lfl2l . 
In their determinant formula for generating function, we can 
set only a single counting parameter, Xk> to be non-zero. 
Then, for arbitrary scattering matrix, the cumulant generat- 
ing function of a current through a corresponding reservoir is 
given by 



S(Xk) = In f l[(C k (E) + A k (E)e Xk + B k (E)e- Xk )j , 

(17) 

where A k , B k , and C k are constants that depend on multi- 
electron scattering amplitudes and channel populations. Ac- 
cording to the golden rule, in the limit of weak transmission, 
coefficients A k and B k can be approximated to first order in 
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s,fc as 



and 



A k (E) =J2ME)(l - n 3 (E))\s kl \ 2 (18) 



B k (E) = - n k {E))n {E)\s ]k \ 2 . (19) 

Here n, (i?) is population of the channel with energy E 1 in the 
i-th contact. Moreover, in this limit, we can use ln(l + x) ~ x 
to approximate 

s(xk) « EE i^i^ 1 - ^>j( e ~ Xk - 1) + 

^^k/n^l-n^e^-l). (20) 

Hence, the counting statistics of individual contact currents of 
non-interacting electrons is equivalent to a sum of two Poisson 
processes, 

S( X k) = A k (e^ - 1) + B k {e-*» - 1). (21) 

This generating function is symmetric under exchange Xk ~> 
~Xk + F k and leads to an FRCC with 



F k = In 



^k \ s kj 



(22) 



C. Networks connected to a single reservoir 

A circuit that is coupled to a single particle reservoir, such 
as shown in Fig.[3ja,b), is another class of systems that show 
FRCCs for currents through some links that, in contrast to pre- 
vious case, belong to loops of the graph. Such networks can be 
used to describe statistics of single molecule events fI51 . For 
example, the graph in Fig. |3ja) corresponds to kinetics of a 
biological enzyme lfl3ll where the external link corresponds to 
the process of enzyme creation/degradation fTD . Probabilities 
of geometric trajectories in such open networks depend only 
on kinetic rates and not on time moments of individual transi- 
tions through links ITD . Consider, for example, the model in 
Fig. [3ja) where we will be interested in currents through the 
link that connects nodes 4 and 1. The arguments in Section 3 A 
lead for this model to a conclusion that it is sufficient to prove 
the FRCC for currents produced by a single particle during its 
life-time because the symmetry of a single particle generating 
function becomes the symmetry of the full counting statistics 
of currents when there is only one external particle reservoir. 
Consider a particle that is just appeared on the node 1 on a 
graph in Fig. [3ja). Let r be the kinetic rate of the transition of 
this particle from node 1 to the reservoir and let fey be kinetic 
rates of transitions between a pair of nodes, i and j. From the 
node 1, the particle can return to the reservoir with probability 
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FIG. 3. (a) and (b): Networks coupled to a single reservoir. Green 
links carry currents that satisfy FRCC and red links carry currents 
that generally do not satisfy FRCC. (c), (d) and (e): Distinct cycles 
of the graph. Cyclic arrows define "+" directions of cycles. 



producing no currents in the system, or it can move into an- 
other node. In the latter case, we know that the particle must 
eventually return to the node 1 because it is the only place 
from which it can escape from the network, and our measure- 
ment time is assumed to be much larger than the particle's 
life-time on the graph. 

Each particle's return to node 1 corresponds either to mak- 
ing none or one of three possible cycles fl3l . shown in 
Fig. |5Jc,d,e). Each cycle can be passed in two directions, 
that we will mark "+" or "— ". Let po, r}\±, r]2±, and 773-1- be 
probabilities of such return events, where indices 1, 2, 3 corre- 
spond to cycles in Fig. [3] respectively, (c), (d), and (e). These 
probabilities depend only on kinetic rates, i.e. they are con- 
stants ifTTl . The generating function of currents through the 
link (4, 1), produced by all geometric trajectories that make 
exactly one return to node 1 is then given by 

ZliXil) = (1 - Pcsc)p C sc[P0 + m + V-3 + 

(m + V- 2 )e Xil + (v-i + m)z~ Xil \- (24) 

After each return to the node 1 the process renews. Hence 
the generating function produced by all trajectories that make 
exactly n returns to the node 1 is given by 

Zn(X4l) = Pesc[(l - Pcsc)(P0 + V3 + V-3 + 

(m + ^ 2 )e X41 + (77-1 + m)e- Xil )T- (25) 

The total single particle generating function is the sum over 
current generating functions, induced by all geometric trajec- 
tories, i.e. 



(26) 



71=0 



which explicitly can be written in a form 

Z(xa) =p esc /[l-(l-p eS c)(AeX"+Be-™+C)}, (27) 

where A = rji + 77_ 2 , B — r)-i + r) 2 , C = p +ry 3 + r/_ 3 . Eq. 
(27 1 is invariant under change of variables, xn ~^ ~Xn + 
F41, where 



F41 = ln([?7_i + rfc]/[r)i + 77-2]), 



(28) 



r/(r + fei2 + fci 6 + feu), 



(23) 



which completes our derivation of the FRCC for the link 
(4, 1). One can easily extend our arguments to all green links 
inFigs.|3j;a,b). 
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IV. EXACTLY SOLVABLE MODELS 

In this section, we derive explicit expression for the gener- 
ating function of currents in the cavity model in several im- 
portant limits and check the presence of FRCC explicitly. We 
show that results are in full agreement with the symmetry de- 
scribed by (|7]l and (111. 



A. Stochastic path integral solution of chaotic cavity model 

To derive FRCCs for the model in Fig. [T] we will em- 
ploy the stochastic path integral technique ll8l ll4|[T5l that was 
previously applied to calculations of current cumulants and 
studies of standard fluctuation theorems in the cavity model. 
This approach is applied to the case of the cavity that has a 
mesoscopic size so that typically we have Q » 1. Follow- 
ing this approach, the counting statistics at a steady state is 
S(x) = Ht, where H = H(x, Qc{x),Xc(x)) is given by 



N 



H = J2 fcj n (Qc)(e*< + ^ - 1) + fc° ut (Qc)(e- 



(29) 

where Q c and \c are expressed through x by solving steady 
state "Hamiltonian equations", 



dH/dQ c = 0, dH/dxc = 0. 



(30) 



To explore those equations, it is convenient to introduce com- 
binations of parameters, 



N 



N 



b (x) = S*exp(- 



(31) 



(32) 



Hamiltonian equations then explicitly lead to relations: 

cxp( Xc ) = [bf 2 /af 1 } 1/2 1 (33) 



/{(^V / ^-«(0))+/2(^V / 77/2-KO))=0, (34) 

where functions /i and f% were defined in Section 2. Gen- 
erally, such nonlinear equations cannot be solved explic- 
itly to determine Q c and x c , but they do imply that Q c , 
as well as combinations ae Xc and be~ Xc , and hence H 
and 5* depend on counting parameters only via the product, 
a(x)b(x)- It is then straightforward to verify that the symme- 
try, S(0, Xk, 0) = S(0, -Xk + Fk, 0, ■ ■ ■), is also 
the symmetry of a(x)b(x) with 



F, = In 



hkY^=£k9t 



k = 1, 



,N, 



(35) 



which proves the FRCC for the cavity model in the limit of 
mesoscopic cavity size. 

The scope of the path integral technique is limited to meso- 
scopic systems. However, we could explicitly verify that the 
FRCC holds true in two exactly solvable limits. In the fol- 
lowing subsections, we show exactly that Eq. ( fTT) ) is satisfied 
for a cavity with exclusion interactions (f 1 = Q max — Q, 
f% = Q)> an d in the classical limit in which the number of 
available states in the cavity is much larger than the number 
of electrons so that in-going rates are not influenced by the 
Pauli principle (/i = 1 and arbitrary ./^(Q))- 



B. Exact solution of cavity model with exclusion interactions 

Here we consider the cavity model with N lead contacts. 
Electrons do not interact with each other except via the ex- 
clusion interactions due to the Pauli principle. Electrons enter 
the cavity through a lead contact i (i — 1, . . . , N) with ki- 
netic rate, k\ n (Q c ) = higi(Q max — Q c ), and leave the cavity 
with rate fc° ut (Q c ) = giQ c . The latter rate is proportional to 
the number of electrons Q c in the cavity. This corresponds 
to non-interacting case. We are interested in statistics of cur- 
rents through a specific lead contact, k. Our exact solution 
of this model is based on the observation that the model is 
equivalent to the model of independent currents through Q max 
quantum dots lfT5lll6l . Each dot can have either zero or max- 
imum one electron inside. If there is no electron in the dot 
then with rates higi, where i — 1, . . . , N, an electron jumps 
into the dot. If the dot has electron inside, then with rates gi, 
i = 1, . . . , N, it leaves to one of the leads. Calculation of the 
cumulant generating function at steady state in such a two- 
state model is straightforward and was discussed in number 
of publications, e.g. fT71 - l20l with a minor difference that here 
we assume many lead contacts. Following e.g. Ifl9l . the cu- 
mulant generating function is given by the larger eigenvalue 
of the matrix 



H 



9k& 



-Xk 



Its eigenvalue, Aq, can be explicitly written as 



(36) 



A = + \/K 2 + 4g k (h k g k + e~^A + eX"h k B), 

(37) 

where A = J2i^k^i9i' B = Ei^feffi' an ^ K is indepen- 
dent of Xk constant. Since the model of the cavity with exclu- 
sion interactions is equivalent to Q max independent processes, 
each having counting statistics Xo(x)t, the cumulant generat- 
ing function for the complete model is given by 



S(Xk) 



Ao(Xk)t- 



(38) 



Obviously, Aq(x) is symmetric under exchange Xk ^ ~Xk 
F k , where F k is given by 
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Fi. = In 



(39) 



in agreement with ( 1 1 



C. Exact solution of the cavity model for stochastic transitions 
with local interactions 



Another exactly solvable model corresponds to a system of 
locally interacting particles performing stochastic transitions 
through a cavity connected to N reservoirs (See Fig. [TJ with 
constant in-going kinetic rates. Let a state vector, \Q C ), be de- 
termined by an occupation number Q c , associated with the 
cavity. Kinetic rates for transitions from the cavity to i-th 
reservoir is given by 



Q c fc(t, Q c )gi, 



(40) 



where f c (Q c ) describes arbitrary local repulsive interactions 
inside the cavity that influence out-going kinetic rates by 
renormalizing single electron free energy in the mean field ap- 
proximation. The in-going kinetic rates through the i-th reser- 
voir is 



which is the result of many -body interactions inside the cavity. 
For the case with local particle interactions, it becomes easier 
to obtain the solutions for the master equation by looking not 
at evolutions of ket-vector but rather at the backward in time 
evolution of bra-vector ( x I' x (t)| given by ifTTl 



dt 



-(* x \C x (t), (* x (t + T)\ = (* x (t)\e- s ^. 

(47) 

According to niHUED, the largest S(t) in (|47]i coincides 



with the cumulant generating function (CGF) of currents per 
period of the parameter evolution in the t — >• oo limit. Follow- 
ing IfTTl . we will search for the solution of ( |47] i in the form of 
a coherent state bra-vector 

(tf x (i)| = (0|exp(C c a c )e- 5(t) , (48) 

and substitute this ansatz in Eq. ( |4"7) , Using the property that 

(* x (*)l4 = <*x(*)lCc, (49) 

and then grouping separately terms near {^ x (t)\f c (Q c )a c and 
separately the remaining functions that multiply (^ x (t)\, we 
find that Eq. ( 47 1 is satisfied if 



(41) 



<c = £^, (50) 
dS/dt = Y,9M t ){ eKi tc-l)- (51) 



Here, for generality, we can also allow arbitrarily prescribed 
explicit periodic time-dependence of ft. = (hi, . . . hjr) 
with period r, keeping other parameters constant. A particle 
distribution function P(Q C ) can be written as a state vector 



* = ^P(Q c )|Qc) 

that satisfies the master equation 

d^/dt = 



(42) 



(43) 



To derive the counting statistics of currents, we should con- 
sider evolution with a twisted master operator, C x that can be 
obtained from the operator C by multiplying its off-diagonal 
elements by factors e ±Xi to count transitions from/to reservoir 
i. For details see Refs.[TT1and[T8l Following those rules, for 
our cavity model, we obtain the form of the twisted master 
operator: 



C x = 9 3 [at{fc{Qc)a c -hje Xl ) 

3 



(44) 



-hj - f c {Qc)a c e x >], Q c = aja c , (45) 



where we have used the "second quantized" version of the 
master equation with 



at|Q c ) = |Q c + l), a c \Q c ) = Q C \Q C - 1) 



(46) 



being the creation/annihilation operators. Here we note that 
operator C x in ( 44 1 is generally non-quadratic in a\ and a c , 



Only parameters hj in |5T| l are time-dependent. It is then triv- 
ial to integrate pT) to find 



S(r, X )=Eft^(e Xi Cc-l), 



(52) 



where we defined hj = h(t')dt'. We can now write 
S(t, x) explicitly as 



S(r, Xk) 



Y, jk hk9k9 3 (^ X3+Xk -^ 



(53) 



To illustrate the presence of the FRCC, we write S(t, x) 
for individual counting fields by setting \j = for j ^ k, 
then 



S(T, X k) = C + B k e-x*+A k e 



Xk 



(54) 



where C = > : ^/' . and 



j = h k 9k Z j9 L k 9j Bk = 9k h9j ^ 



Ej9j 



Ej9j 



The form of CGF in Eq. ( 54 1 has the symmetry under the ex- 
change of Xk -Xk + F k, where F k is given by 



F k = In 



T,j^ k 9j h j 
hk Ef^fcffj, 



(56) 



This expression for F k is the same as ( 1 1 
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V. NUMERICAL CHECK OF FRCC FOR NETWORKS 



ing set of linear equations: 



In this section, we consider more complex networks and 
analyze the validity of the FRCC for any particular link nu- 
merically. 



A. Network with loops and backbone link 



Cfe) - 9k((k 



~Xk 



) = 0. (59) 



S( Xk )-S(- X k+F k ) = 




Lbfc 1-2 




Link 1-3 




Link 1-5 




Link 3-5 




Link 2^ 




Link 2-6 




Link 4- 6 





Explicit functional form of S, even for a circuit of a moder- 
ate size in Fig. [2] is quite bulky to be shown here. In any case, 
it is difficult to observe the presence of an FRCC just by look- 
ing at an analytical expression for a generating function. To 
demonstrate the presence of an FRCC, we resort to numerical 
solution of ( 

m 



in 



We set all but one counting parameters 
to be zero. The presence of an FRCC can 
be checked by looking at the solution contours of the equation 
S(x)-S(-x+F) = plotted as a function, F(x). An FRCC 
occurs for the link k if there exists a solution contour for which 
F k is independent of Xk- In other words, if we plot all pairs, 
(Xk,F k ), that satisfy equation, S(Xk) ~ S(~Xk + F k ) = 0, in 
plane with axes Xk and Fk, then if an FRCC holds, the curve 
must be a horizontal line, which is parallel to x^-axis. 

The plot in Fig. [4] shows a solution Fij(x%j), where 
runs through all internal links of the network in Fig. [2] It 
clearly shows that only for the link (1,2) we have F12 = 
const. For all other links, solution contours are nontrivial 
functions (-Fy (x-h) ^ const). The FRCC is upheld by the 
links that do not belong to any loop. Such a link in our ex- 
ample is the (1, 2) link, and this is the only link that supports 
FRCC, which is in agreement with our discussion in the main 
text. 



FIG. 4. Numerical plots for the solution contours of S(x) — S(~x + 
F) — for currents in geometry of Fig. [2] Parameters were set to 
numerical values: gi2 = 0.3854, 313 = 0.6631, gin = 0.5112, 
ff 3 5 = 0.62 53, #24 = 0.17483, g 26 = 0.02 5 97, g 46 = 0.93 1946, 
Ski = 9ii, 9i = 0.5379, g 2 = 0.004217, g 3 = 0.8144, g 4 = 
0.8396, g 5 = 0.1765, g 6 = 0.8137, hi = 8.7139, hi = 8.8668, 
h 3 = 4.6144, h 4 = 8.1632, h 5 = 7.6564, and h 6 = 8.01407. 

Consider a graph that is coupled to external particle reser- 
voirs and has loops and a backbone link as shown in Fig. [2] 
For numerical check we restrict to models with only one reser- 
voir per node. Generalizing the cavity model we define addi- 
tional parameters that are the kinetic rates of transitions from 
node i to node j, given by fey = Qifigij. The twisted master 
operator for this case can be written as ifTTI 



At = ^9jka]{e Xlk f k a k - fjUj) 



jk 



(57) 



where we introduced additional parameters Xij t° count cur- 
rents through internal links Following the prescription 
for the case of the single cavity (node) model we obtain the 
CGF, S(xjk, Xk), which reads: 



s(t,x) = X)fiM^-i). 



(58) 



where parameters £'-s can be obtained by the solving follow - 



B. Numerical check for networks coupled to a single reservoir 

In this subsection, we numerically analyze the class of net- 
works with loops. Particles can enter and leave only via a 
single reservoir. Figs. |3ja,b) show networks, in which all in- 
ternal links belong to some loops of a graph. Detailed balance 
on kinetic rates is not assumed. The twisted master operator 
for this case is 



^9]ka){e X]k f k a k - fcaj) 



jk 

- 9n[a ] n {f n a n - h n e Xn ) +h n - f n a n e~ Xn ]. (60) 
To account for the lack of detailed balance condition we 



simply allow giy 7^ g^. In Eq. (60 1 the index "n" corresponds 



to the only node in the cyclic network that is connected to 
the reservoir. Similarly to the previous example, we obtain 
the following cumulant generating function for the case of the 
networks in Figs. [3ja,b), 



S(r, x) = gnh„(e Xn C n - 1), 



(61) 



where the average occupation number of the node Q can be 
obtained by solving the following set of linear equations, 



Ck) — <5/c,n5fe(Cfc 



-Xk 



; ) = 0. (62) 



The numerical check for this case is done by inspecting the 
solution contours of S( X ) — S(— X + F) = 0, as described in 
previous subsection. 
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S(xk)-S(- Xk +F k ) = 



2 



4 







Li* 1-2 = 

Lint 3-3 

Liilcj-J 

- . i - 










Lille [-4 
Link 5-6 










Link 6-1 ^— 








10 -5 ( 


1 5 10 



Xk 



FIG. 5. Numerical plots for the solution contours of S(x) — S(—x + 
F) — for all internal links in Fig. |3ja). We consider parameters 
having the following values: 312 = 0.4557, 321 = 0.6282, 323 = 
0.97782, 332 = 0.04 8 59, 334 = 0.9432, 343 = 0.7787, 345 = 
0.9622, g 54 = 0.42 98, g 56 = 0.30 23, g 65 = 0.38 56, fifei = 0.4667, 
gie = 0.06163, 341 = 0.2779, 314 = 0.09021, 31 = 4.7019, 
hi = 8.7658. 




Li* 1-2 




Unit 2-3 




Link 1-1 




Lmk-l-S 




Link 1-4 




Link 5-6 




Link 6-1 





FIG. 6. Numerical plots for the solution contours of S(x) — S(—x + 
F) = for all internal links in Fig. [5Jb). The network parameters 
for this case are: 512 = 4.5654, 521 = 7.83, (723 = 8.6823, 332 = 
0.6629, 334 = 7.0427, 343 = 7.4394, 345 = 7.95, 354 = 0.1431, 
g 5 6 = 0.40 52, yes = 2.5126, g 61 = 9.5782, g 16 = 0.08372 ,341 = 
9.0737, 3i4 = 0.81007, 32 = 5.5299, and hi = 8.2477. 



We consider two different cases: 

(a) when the reservoir is connected to the central node 1, 
and 

(b) when the reservoir is connected to a non-central node 
(e.g. the node 2). 

We first consider case (a). According to the plot (see Fig.BJ, 
all the links have solution contours that are independent of x- 
Hence the FRCC holds for all the links when the reservoir 
is attached to the three link junction 1. The degeneracies in 



the values of F are due to the charge conservation. Currents 
through the links (1, 2), (2, 3) and (3, 4) are the same which 
leads to the same value of F. This degeneracy is also seen for 
the links (4, 5), (5, 6) and (6, 1). Hence, in Fig. [5] we see only 
three sets of degenerate horizontal lines. 

We perform similar numerical analysis for the case (b) 
when the reservoir is connected to node 2. The solution con- 
tours for S(x) — S(~x + F) = are shown in Fig. [6] For this 
case we obtain degenerate constant lines of F — const for the 
links (1, 2), (2, 3) and (3, 4), thereby satisfying an FRCC. All 
the other links do not satisfy an FRCC. 



VI. CONCLUSION 

We demonstrated that FRCCs appear in fundamental mod- 
els of nanoscale electric circuits. Many-body electron inter- 
actions, including Coulomb interactions and exclusion inter- 
actions, do not break the FRCC prediction in the model of a 
chaotic cavity coupled to several leads. This robustness can be 
used to extract information about relative sizes of single parti- 
cle tunnelling barriers independently of electron interactions. 
To the best of our knowledge this is a unique example in which 
measurements of fluctuation relations can provide quantitative 
information that cannot be easily obtained by measuring low- 
est current cumulants at given voltages. We also demonstrated 
that the FRCCs extend to the quantum regime of coherence 
among electron trajectories. FRCCs can be exact even when 
generating functions cannot be derived. This reflects the fact 
that FRCCs follow from the properties of single particle geo- 
metric trajectories that separately show profound symmetries 
even when the complete stochastic evolution of a system is 
very complex. 

An unusual property of FRCCs is that they are not directly 
related to the system's dissipation function. This distinguishes 
them from the vast number of previously found fluctuation re- 
lations for currents, entropy, heat and work. Instead, the fact 
that FRCCs follow from constraints on geometric trajectories 
relates them to the principle of Geometric Universality of Cur- 
rents 1 1 1 1 and the class of exact results in non-equilibrium 
physics called "no-pumping theorems" lfTTll2~Tl . On the other 
hand, similarity between FRCCs and standard fluctuation re- 
lations suggests that there can be a more fundamental theory 
in the background of both fluctuation and no-pumping rela- 
tions. 

Our results should stimulate further theoretical and experi- 
mental studies of fluctuation relations, including the search for 
unifying fundamental principles in non-equilibrium statistical 
mechanics and new measurement techniques that are enabled 
by FRCCs. 

Recently, we have learned of another work that obtained a 
large class of FRCCs ll22l that do not directly relate to the 
dissipation functional of the total system. Unlike our results, 
their FRCCs apply to a closed network topology without ex- 
ternal particle reservoirs. 

Acknowledgment. Authors thank D. Andrieux for useful 
discussion. The work at LANL was carried out under DOE 
Contract No. DE-AC52-06NA25396. Work at NMC was sup- 



10 



ported by NSF under ECCS-0925618. 



[1] G. N. Bochkov and Yu. E. Kuzovlev, Sov. Phys. JETP 49 543 
(1979); Physica A 106, 443 (1981); D. J. Evans, E. G. D. Co- 
hen, G. P. Morriss, Phys. Rev. Lett. 71, 2401 (1993); C. Jarzyn- 
ski, Phys. Rev. Lett. 78, 2690 (1997) 

[2] M. Esposito, U. Harbola, and S. Mukamel, Rev. Mod. Phys. 81, 
1665 (2009); P. Hiinggi and F. Marchesoni, ibid. 81, 387 (2009) 

[3] D. Andrieux and P. Gaspard, J. Stat. Phys. 127, 107 (2007); H. 
Forster and M. Buttiker Phys. Rev. Lett. 101, 136805 (2008); 
D. Andrieux, P. Gaspard, J. Stat. Mech. P02057 (2009); N. A. 
Sinitsyn, A. Akimov, and V. Y. Chernyak, Phys. Rev. E 83, 
021107 (2011); J. Ren, V. Y. Chernyak and N. A. Sinitsyn, J. 
Stat. Mech. P05011 (2011); R. Sanchez, R. Lopez, D. Sanchez, 
and M. Buttiker, Phys. Rev. Lett. 104, 076801 (2010); A Rakos 
and R J Harris J. Stat. Mech. (2008) P05005 

[4] M. Esposito, U. Harbola, and S. Mukamel, Phys. Rev. B 75, 
155316 (2007); Y. V. Nazarov, Y. M. Blanter, "Quantum Trans- 
port", Cambridge University Press, New York (2009) 

[5] Gustavsson et al, Phys. Rev. Lett. 96, 076605 (2006); Sukho- 
rukov et al, Nature Phys. 3, 243 (2007) 

[6] Y. Utsumi et al, Phys. Rev. B 81, 125331 (2010) 

[7] G. B. Cuetara, M. Esposito, P. Gaspard, Preprint 
arXiv/1 105.5974 (2011) 

[8] S. Pilgram et al, Phys. Rev. Lett. 90, 206801 (2003); A. Al- 
tland, A. De Martino, R. Egger, and B. Narozhny, Phys. Rev. 
Lett. 105, 170601 (2010) 

[9] P. W. Brouwer, Phys. Rev. B 63, 121303 (2001) 



[10] P.-E. Roche, B. Derrida and B. Douot, Eur. Phys. J. B 43, 
529541 (2005); A. N. Jordan, E. V. Sukhorukov, and S. Pil- 
gram, J. Math. Phys. 45, 4386 (2004) 

[11] V. Y. Chernyak, M. Chertkov and N. A. Sinitsyn, J. Stat. Mech. 
P09006 (2011) 

[12] L. S. Levitov, and G. B. Lesovik, JETP Lett. 58, 230 (1993) 

[13] D. R. Astumian, Annu. Rev. Biophys. 40 289313 (2011) 

[14] N. A. Sinitsyn, N. Hengartner, and I. Nemenman, PNAS, 

U.S.A. 106, 10546 (2009) 
[15] N. A. Sinitsyn and I. Nemenman, Phys. Rev. Lett. 99, 220408 

(2007) 

[16] J. Ohkubo, Phys. Rev. E 80, 012101 (2009) 
[17] B. Derrida and J. L. Lebowitz, Phys. Rev. Lett. 80, 209 (1998) 
[18] D. A. Bagrets and Y. V. Nazarov, Phys. Rev. B 67, 085316 
(2003) 

[19] I.V. Gopich and A. Szabo, J. Chem. Phys. 124, 154712 (2006) 
[20] N. A. Sinitsyn and I. Nemenman, Europhys. Lett. 77 58001 
(2007) 

[21] D. Astumian, Proceed. Nat. Acad. Sci. U.S.A. 104, 19715 
(2007); S. Rahav, J. Horowitz and C. Jarzynski, Phys. Rev. 
Lett. 101, 140602 (2008); J. E. Horowitz and C. Jarzynski, J. 
Stat. Phys. 136, 917 (2009); V. Y. Chernyak and N. A. Sinitsyn, 
Phys. Rev. Lett. 101 160601 (2008); C. Maes, K. Netocny, S. 
R. Thomas, J. Chem. Phys. 132, 234116 (2010). 

[22] A. C. Barato, R. Chetrite, H. Hinrichsen, and D. Mukamel, 
Preprint: arXiv/1 109.2517 (2011) 



